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Abstract A novel approach to electronic correlations and magnetism of crystals 
based on realistic electronic structure calculations is reviewed. In its 
simplest form it is a combination of the "local density approximation" 
(LDA) and the dynamical mean field theory (DMFT) approaches. Us- 
ing numerically exact QMC solution to the effective DMFT multi-orbital 
quantum-impurity problem, a successful description of electronic struc- 
ture and finite temperature magnetism of transition metals has been 
achieved. We discuss a simplified perturbation LDA+DMFT scheme 
which combines the T-matrix and fluctuation-exchange approximation 
(TM-FLEX). We end with a discussion of cluster generalization of the 
non-local DMFT scheme and its applications to the magnetism and su- 
perconductivity of high-T c superconductors. 
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1. Introduction 

The theory of electronic structure and magnetism of solids histori- 
cally was split into two distinct parts, namely, the model investigations 
of many-body effects and the calculations of the energy spectra and prop- 
erties of specific compounds in the framework of density functional (DF) 
scheme [1,2]. Recently, within the dynamical mean-field theory (DMFT, 
for a review see Ref. [3]) the correlation effects have been incorporated 
into realistic electronic structure calculations [4-10]. This method has 
been successfully applied to a number of classical problems of solid state 
physics such as the finite-temperature magnetism of iron-group metals 
[11], a — 5 transition in plutonium [12], electronic structure of doped 
Mott insulators [9]. In contrast with standard DF theory, in this new 
approach known as "LDA+DMFT" [4] or "LDA++" [5] the total en- 
ergy of the system (or, more accurately, the thermodynamic potential 
f2) is considered as a functional of the Green function instead of the 
denstity matrix [7, 8, 13, 14]. To stress this new feature more explicitely 
we will use the term "spectral density functional (SDF)". The analyti- 
cal properties of the Green function garantee that the knowlegde of the 
spectral density is equivalent to the knowledge of the time-dependent 
Green function whereas the density matrix is just static value of the 
latter [15]. Here we will describe the basic ideas of the SDF method, 
both in the framework of a standard DMFT and from a more general 
point of view, discuss possible cluster generalizations of the DMFT, and 
consider the applications of SDF to correlation effects and magnetism in 
transition metals. 

2. Dynamical Mean Field Theory: an effective 
action perspective 

A most economical approach to unifying the various dynamical mean- 
field approximations in use, is provided by the effective action construc- 
tion [16]. The idea is to select a set of variables which is relevant to the 
physics of the problem and to write down a functional of the relevant 
variables. The extremum of this functional yields the values of those 
variables in equilibrium, and the value of this functional at stationarity 
gives the free energy of the system in equilibrium. 

Density functional theory is the simplest example of this construction 
[17], here the total energy of the solid is expressed in terms of the density 
of the electrons. Another well tested example is the spin density func- 
tional theory in which the total energy is expressed in terms of the spin 
and charge densities. The construction of an explicit expression of the 
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exact effective-action functional is usually not available except for the 
case of very simple examples and the success of the method relies on the 
availability of good approximations to this functional. The LDA and 
the LSDA approaches have been extraordinarily successful for weakly 
correlated systems. 

The dynamical mean-field approach to model Hamiltonians of strongly 
correlated systems on a lattice can also be brought to such an effective 
action perspective. More important this construction can be easily gen- 
eralized to incorporate the so-called extended dynamical mean-field ap- 
proach (E-DMFT) [18], by constructing an effective action for both the 
local Green function and the local density-density and local spin-spin 
autocorrelation function. The effective-action approach also allows the 
formulation of the E-DMFT approach in the case where spatial and spin 
or charge symmetries are broken. [8, 19]. 

The effective-action approach, allows us a simple combination of the 
density functional theory and the DMFT. By Legendre transformation 
techniques one can construct a functional of the density and of the local 
spectral function of the heavy orbitals whose extremization would yield 
the exact density local spectral function of the heavy orbitals, and the 
total energy of the system. Again the exact form of this functional are 
not known in explicit form, but useful approximations to it are available, 
and can be used to study interesting problems. 

This realistic DMFT, LDA+DMFT approach or LDA++ approach, 
was first implemented ignoring the coupling between the density and the 
local spectral function. This amounts to performing first a LDA calcula- 
tion to derive a tight-binding model Hamiltonian, and then performing 
a DMFT calculation for the spectra of the DMFT Hamiltonian. This is 
very close in spirit to the philosphy of model Hamiltonian calculations. 
Recently a full implementation of the self-consistent determination of 
the density and the local spectra was carried out by S. Savrasov [12], 
and this is now closer in spirit to traditional first-principles electronic 
structure calculations. 

The effective action approach, can be generalized to clusters, if short- 
range correlations need to be taken into account. In the context of model 
Hamiltonians, a functional of the restriction of the Green function to a 
given cluster is defined. The extremization of this functional give rise 
to cluster dynamical mean field equations. There is no difficulty in con- 
structing hybrid funtionals for the self-consistent determination of both 
the density and the cluster Green functions, namely a CDMFT+LDA 
method. 

In the next section we motivate the approximate DMFT form of those 
functionals from a perspective of a reduced fermionic description. 
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3. Fermionic reduced description and dynamical 
mean-field theory 

In the SDF approach two-electron subsytems are introduced, one of 
them is described by the standard DF theory (usually this is the subsys- 
tem of sp-electrons) and dynamical interelectron correlations are taken 
into account for the another one (usually this is the subsystem of d- or /- 
electrons). Further simplifications can be connected with the local form 
of interelectron interactions (only Hubbard-type on-site electron corre- 
lations are considered), using some approximations like DMFT or its 
cluster generalizations [22, 23, 24], etc. From a general point of view all 
these approaches can be considered as specific cases of a "coarse grain- 
ing" (reduced description) ideology [25] when all the variables describing 
the system can be separated into the "gross" and "slave" variables; the 
only assumption that a closed set of the equations of motion for the gross 
variables exist is sufficient to find an explicit form of these equations [26] . 
Here we demonstrate the coarse graining procedure for fermionic degrees 
of freedom for the system of interacting electrons in a crystal. 

Let us start with the functional integral over the Grassman creation 
and annihilation electron fields c + ,c, where the "measure" is exp (— S), 



+ S mt [c + ,c] (1) 



Here n is the chemical potential, (3 = 1/T is the inverse temperature, and 
Hq is the "one-particle" part of the Hamiltonian. In the SDF approach 

H = - V 2 /2 + V KS (2) 

where Vks is the Kohn-Sham self-consistent potential [1]. 

In the spirit of a reduced description approach [25] we introduce 

"gross" variables d = f • c (more explicitly, d^ = f^j c j where j are site 

j 

indices; e.g., in a standard DMFT f^j = S^Sjo). The effective action 
for the gross variables S e ff is defined by introducing a delta-functional 
6(d — f ■ c)8 (d + — f* ■ c + ) [27] into the functional integral. 



S 



S ix S 



dr 



„+ 
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exp (S e ff) = — J DcDc + exp(-S) * 

J D\D\ + exp [iX (d + - f* ■ c+) + iX + (d - f ■ c)] (3) 

Unfortunately, this functional integral cannot be estimated exactly 
and approximations are needed. Here we consider only the simplest ap- 
proximation to the derivation of the effective action, we replace Sj n t [c + , c] 
Sj^[d + ,d] but compensate for the omission of the interaction terms 
away from those in the gross variables, namely in the medium by adding 
a medium self-energy in the standard DMFT 

S^S = S + Y,f I dTdT'c + (T,i)Z(T-T',i,i)c(T',i) (4) 

i>0 







In other words, the DMFT approach treats exactly the interaction 
terms only in those part of the interaction Lagrangian which can be 
written in terms of gross variables, while the rest of the interaction terms 
are handled in a gaussian approximation by replacing the interaction 
terms by a medium self-energy. 

Now the effective measure for the gross variables is given by 

exp (-S eff ) = — J DcDc + exp (-S C G + * 

DXDX+ exp [iX (d + - f* ■ c+) + iX + (d - f ■ c)] (5) 



/ 



where Z = J DcDc+ exp (-Sea) , S C G = S med + Sgf [d+ , d] . 
Passing to the Matsubara frequencies [15] one has 

Smed = -T^2c + (iu n j) [iu n + [I - H - S (iu n )\ c {iuJnj) (6) 

In the standard DMFT scheme the correlations between d(/)-electrons 
on a central site j = are treated exactly while the correlations in the 
medium are treated in the Gaussian approximation. It is sometimes 
stated in the literature that this assumption violates the translational 
invariance of the system. The effective-action approach, is meant to 
construct a functional of the selected variables, so by construction it is 
not possible to discuss within this framework the issue of translation 
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invariance. The "self-consistency condition" does the best possible job 
within a Gaussian approximation to restore the equivalence of the central 
site and the medium. 

We show below that this reasoning in a cluster setting leads to the 
"cellular DMFT" approach [24]. In principle, the correlations between 
the gross and slave variables (e.g., the long-range part of the Coulomb 
interaction) also can be taken into account in a close analogy with the 
classical spin models [25]; the corresponding modification of the DMFT 
approach will be considered elsewhere. In accordance with the general 
scheme of the non-equilibriun statistical operator method [26], we add 
an "auxiliary field" conjugated to the Green function of the gross vari- 
ables, namely, the term d + T,d. These fields will restore the translational 
invariance for the fermionic Green function (see below, Eq.(ll)). 

Calculating the Gaussian integral over c, c + one has 

exp (S eff ) = exp (-Sgf [d + , d] + d + td) * 

J D\D\ + exp [i\d + + i\ + d - X + f ■ G med ■ f*X] (7) 

where the medium Green function reads 

G med = [iuj n + n - H - S {iuJnYT 1 ■ (8) 
At last, calculating the Gaussian integral over A,A + we obtain 

Seff [d + , d] = Sgf [d + , d]-TY,d + {iu n ) So" 1 (uj n ) d (iu n ) (9) 

where 

Go 1 (»wn) = (/ ■ G med ■ f* r 1 + E (10) 

The self-consistency condition that determines S requires the equality 
of the local Green function computed from the reduced description and 
from the medium agree, i.e., 

(Td(r)d + (r')) s = (Td(r)d + (r')) Seff (11) 

Further, (Tc(t)c + (t')) s = G me d (because of the replacement Si n t [c + ,c] - 
Sgf [d+,d])and 

(Td(T)d + (T f )) s = f-G med -f* (12) 
On the other hand, (Td(r)d + (r')) s = Q satisfies the Dyson equation 



<r 1 = g, 1 - s CG = (/ • c med ■ ry 1 + s - s CG (13) 
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and we have S = £cg where Scg is a "coarse-grained" self-energy. The 
most natural choise is T*cg = /*■£■/. 

To obtain the main equation of the DMFT [3] we have to choise / as 
a projection operator on the central site. In this case Eq.(13) will take 
the desired form 

-l 



G 1 (iu n ) = 



E 



+ £ (iLO n ) 



(14) 



k iw„ + /i - i (k) - S (i(jj n ) 

where the quasimomentum k runs the Brillouine zone, t (k) is the Fourier 
transform of the Hamiltonian Ho projected into the subspace of d(f)- 
electrons. 

To consider posiible cluster generalizations one can choose let / = 
PU + P where P is the projection operator on the cluster and U is an 
unitary transformation of the variables in this cluster. Thus our result 
reads: 



Go 1 



iu n + fj, - t (kc) - S 



+ c/S?7 + (15) 



where k c runs the "new" supercell Brillouin zone, L and L c are the sizes 
of the crystal and cluster, correspondingly. In terms of Scg = UT>U + , 

-l 

| \ - ' 



G 1 (i^n) 



r- ) |^ia; n + ^-C/t(k c )?7+-S CG 



+ (16) 



As it was shown in Ref. [24] this equation for general basis set can be 
very useful in the optimisation of interaction problem within the cluster 
or cellular DMFT scheme. Another version of the cluster DMFT [23] 
will be described below. 



4. Spectral density versus density functionals 

In a standard DF theory the thermodynamic potential for non-correlated 
conduction- "c" electrons £l c is represented as a functional of the electron 
density p (r) which is, generally speaking, a matrix in spin indices. For- 
mally it can be represented as a thermodynamic potential of the Kohn- 
Sham quasiparticles [1], 0, sp , minus the contribution of the so called 
"double counted" terms, Qdc'- 



O c — O c 



-Tr \og[iuj + V 2 /2 - Vks] (17) 
J V KS (r)p(r)dv- J V ext (r)p(r)dr-±J P ^ M ^ dvdv' - E xc [p] 
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where Tr = Tr^iLa , Tr^ is the sum over Matsubara frequencies Tr w . . . = 
TY,e luj0+ uj = 7rT(2n + 1) , n = 0,±1, T is the temperature, and 

iLa are site numbers (i), orbital quantum numbers (L = l,m) and spin 
projections a, respectively, V ex t(r) is the external potential, E xc [p] is the 
exchange-correlation energy, and the Kohn-Sham effective potential is 
defined as 

^( r ) = W r) + /^ r < + f^. (18, 

In contrast with the standard density functional theory, the SDF ap- 
proach deals with the real dynamical quasiparticles for correlated "d- 
electrons" defined via local Green functions rather than with Kohn-Sham 
"quasiparticles" which are, strictly speaking, only auxiliary states to cal- 
culate the total energy. Therefore, instead of working with the thermo- 
dynamic potential 0, as a density functional we have to start from its 
general expression in terms of an exact Green function [28] 



Q d sp = -Tr {in [S-Gq- 1 ]} 

nj c = TrSG-$ (19) 

where G, Go and S are an exact Green function, its bare value and 
self-energy, $ is the Luttinger generating functional (sum of the all con- 
nected skeleton diagrams without free legs), respectively. A complete 
SDF thermodynamic potential is equal to Q = tl c + Q d . We have to keep 
in mind also the Dyson equation 

G- 1 = Gq 1 - S (20) 

and the variational identity 

s =«r < 21 > 

When neglecting the quasiparticle damping, £l sp will be nothing but the 
thermodynamic potential of "free" fermions but with exact quasiparticle 
energies. Formal analogies between Eqs.(17) and (19), (18) and (21) 
are obvious: the self-energy plays the role of the Kohn-Sham potential 
(without the external potential) and the Green function plays the role of 
the density matrix. As an example of this correspondence one can prove 
[7] in the framework of the SDF an useful identity known as the "local 
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force theorem" basically in the same way as it has been done within DF 
theory [29, 30]. 

For both parts of the thermodynamic potential, £l c and Q d , the local 
force theorem is based on extremum properties with respect to the vari- 
ation of the density matrix and Green function, respectively. Further we 
will consider the contribution Q d (omitting the index d for brevity). In 
principle, fermionic reduced description scheme allows us to combined 
conduction "c" and correlated "d" part of the total SDF: 9. = tt c + tt d 
[14]. 

5. Effective exchange interactions 

Let us discuss the problem of calculation of effective exchange in- 
teractions (Jij) in correlated systems. In principle the Jy parameters 
are not well defined for arbitrary magnetic systems, and the traditional 
way to study spin excitations related to the calculation of non-local 
frequency-dependent spin-susceptibility [3, 32]. In this case the auxiliary 
space-time dependent magnetic field is added to the the Hamiltonian: 
<rh(r, r) and the second derivative of the free-energy with respect to the 
magnetic field gives the interacting spin-susceptibility: x" 1 = Xq 1 ~ T, 
where xo is an empty-loop susceptibility and T is the vertex correction 
[3, 32]. Here we consider a simple approximation of "rigid spin rotation" 
of spectral density for a small angle: 

5ei = Sifi x e« (22) 

where is a general direction of constrained effective spin-dependent 
potential on site i and 8<pi is a rotation vector. In this case it is use- 
ful to write explicitly the spinor structure of the self-energy and Green 
functions: 

= T.i + ^a (23) 
Gij = G\j + Gfj<j 

where T,\ c,s ^ = \ ± S| = Sfe^, a = (a x ,a y ,a z ) are Pauli ma- 
trices, G°j = ^Tr a (Gij) and Gf^ = ^Tr a (Gija). We suppose that the 
bare Green function G° does not depend on spin direction and all the 
spin-dependent terms including the Hartree-Fock terms are incorporated 
in the self-energy. In the rigid spin approximation we assume that the 
unit vector ej does not depend on the energy and orbital indices and 
represents the direction of the average local magnetic moment on the 
site i. Note, that the thermodynamic potential O should be considered 
as a constrained SDF which depends on as on external parameters 
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(cf. Ref. [33]). Then the variation of the thermodynamic potential with 
respect to small spin-rotation can be written as 

5£l = 5*Q sp + 5i£l sp — 5£ldc (24) 

where 5* is the variation without taking into account the change of the 
"self-consistent potential" (i.e. self-energy) and S\ is the variation due 
to this change of S. Taking into account Eq. (21) it can be easily shown 
(cf. Ref. [28]) that 

S&sp = dn dc = TrGSY, (25) 

and hence 

(26) 



SSI = 6*n sp = -S*Tr In 



£ - Gq 1 



which is an analog of the "local force theorem" in the density functional 
theory [30]. 

In the case of rigid spin rotation the corresponding variation of the 
thermodynamic potential can be written as 

6(1 = Vrftpi (27) 

where the torque Vj is equal to 

V, = 2T ru)L [Ef x GfJ (28) 

Based on the expansion of this expression (28) in a sum of pairwise 
contributions one can obtain [7] useful formula for the effective magnetic 
interactions: 

Jij = -Tt^l (S?4S|4) (29) 
and, correspondingly, for the stiffness tensor of a ferromagnet: 

^--->^{^^) <-> 

where M is the magnetic moment per unit cell. These results generalize 
the LSDA expressions of Ref. [30] to the case of correlated systems. 
Note that passing from Eq.(28) to Eq.(29) is not accurate, since the 
exchange parameters are connected with the second variations of the f2- 
potential and use of the local force theorem can not be justified. Eq. (29) 
corresponds to the "empty loop" approximation neglecting the vertex 
corrections. At the same time, for the stiffness the latter are absent and 
Eq.(30) appears to be exact provided that the self-energy and three-leg 
vertex are local (as in the DMFT) [13]. 

The fact that vertex corrections to the spin stiffness are absent within 
DMFT, is also suggested by the analogy between electric and spin trans- 
port developed in Ref. [20]. The spin- wave stiffness can be obtained from 
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the zero-frequency limit of the spin conductivity. Within DMFT , the 
charge current-current correlation function does not require vertex cor- 
rections [21] and can therefore be obained directly from the convolution 
of the one electron spectra. These arguments are independent of the 
spin structure, and can therefore be used for the spin-wave stiffness. 

In order to elucidate the approximation behind the expression for 
the exchange parameters (Eq. (29)), we consider the energy of a spi- 
ral magnetic configuration with the rigid rotation of the spinor-electron 
operators by the polar angles 9 and (p: 

C-im > U (6i, Ifii) Ci m 

where 

cos6>/2 sin 9/2 exp (—if) \ 

-sin 0/2 exp (i^) cos 9/2 ) ^ ' 

assuming that 0j = const and tpi = qRi where Ri is the site-lattice vector. 
Since we take into account only on-site correlation effects the interaction 
term in the Hamiltonian is invariant under that transformation, and the 
change of the Hamiltonian is 



u(e,<p) 



6H = J2Tr n 



1 c, 



S ± H + 5 2 H 



SxH = sin 2 |E Tr ma (t (k+q) - t (k)) c+c k 

k 



(32) 



S 2 H= ± sin #£Trv f 



(exp (iqRj) - exp (iqRj)) 



Consider further the case of small 9, we can calculate change of the 
total energy to lowest order in 6 corresponds to the first order in 5\H 
and the second order in 5 2 H: 



SE = j{E[t(k+q)-t(kfrffr m |^j[<(k+q)-t(k)] 
7 (k, q) G l (k + q) [t (k+q) - t (k)]G T (k)}, 



where = Tr ma (^c^Ck/, q is a four- vector with component (q,0), and 
7 is the three-leg vertex. Our main approximation is to neglect of the 
vertex corrections (7 = 1). In this case the previous equation takes the 
following form: 



SE = -^Tr m {.|^[t(k+q)-t(k)]* (33) 
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G i (k+q) [GJ 1 (k + q)- Gf 1 (k) + i (k+q) - t (k)] G T (fc) } 
Using the following consequence of the Dyson equation: 



t (k+q) - t (k) = Gf 1 (*) - G7 1 (fc + ?) + S T (E) - E 4 (£) (34) 

one can rewrite Eq. (33) in the form: <5i£ = ^-[J(O) — ■/(<?)] with the 
exchange integrals corresponding to Eq. (29). We conclude that the 
expression for Jij is accurate if the vertex corrections can be neglected. 
Note that in the limit of small q this can be justified rigorously, provided 
that the self-energy and three-leg scalar vertex are local. Therefore, 
the expression for the stiffness constant of the ferromagnet (Eq. (30)) 
appears to be exact in the framework of DMFT [13]. 

It should be stressed that the exchange integrals discussed here are 
just static characteristics connected with energy of inhomogeneous spin 
configurations. They determine the frequencies of spin excitation in 
itinerant-electron systems only under adiabatic (rigid-spin) approxima- 
tion. For more general consideration one should calculate q- and in- 
dependent spin susceptibility. 

6. Electron correlations and finite-temperature 
magnetism in transition metals 

Now we describe the applications of the SDF approach to a classi- 
cal problem of finite-temperature magnetism of the iron-group transi- 
tion metals. Despite a lot of attempts starting from seminal works by 
Heisenberg and Frenkel (for review of early theories see e.g. [31, 34, 35]) 
we have yet no complete quantitative theory describing their magnetic 
and spectral properties. The reason is that to describe the properties of 
Fe, Co, and Ni one has to solve the problem of taking into account mod- 
erately strong electronic correlations where approaches developed both 
for weakly correlated systems such as normal-group metals and to higly 
correlated systems such as Mott insulators are not, generally speaking, 
reliable. There were many attempts to introduce by some way correlation 
effects in band structure calculations of these metals [36, 37, 38, 39, 40], 
not referring to numerous purely model works. But the question of 
applicability of specific approximations such as lowest order perturba- 
tion theory [37, 40], moment method [39] or three-body Faddeev equa- 
tions [38] is not clear and one needs a reliable approach which would be 
checked carefully for model systems and demonstrate its applicability for 
moderately correlated systems. It has been demonstrated in Ref. [11] 
that the ab initio dynamical mean-field theory does give a very success- 
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ful description of both correlation effects in the electron energy spectra 
and the finite-temperature magnetic properties of Fe and, especially, Ni. 
Here we present the corresponding results. 

We start with the LDA Hamiltonian in the tight-binding orthogonal 
LMTO representation H^^(k) [41], where m describes the orbital basis 
set containing 3d-, 4s- and Ap- states, and k runs over the Brillouin zone 
(BZ). The interactions are parameterized by a matrix of screened local 
Coulomb interactions, U mm i, and a matrix of exchange constants, J mTO ', 
which are expressed in terms of two screened Hubbard parameters, U 
and J, describing the average Coulomb repulsion and the interatomic 
ferromagnetic exchange, respectively. We use the values U = 2.3 (3.0) 
eV for Fe (Ni) and the same value of the interatomic exchange, J = 0.9 
eV for both Fe and Ni, a result of constrained LDA calculations [42, 4, 5]. 
These parameters, which are consistent with those of many studies result 
in a very good description of the physical properties of Fe and Ni. 

As was discussed above, the DMFT maps the many-body system onto 
a multi-orbital quantum impurity, i.e. a set of local degrees of freedom in 
a bath described by the Weiss field function Q. The impurity action (here 
n rn u = cf na c m(7 and c(r) = [c mo -(r)] is a vector of Grassman variables) 
is given by: 

S eff = - [ dr f dT'Tr[c + {r)g-\T,T')c{T')\ + 
JO Jo 

\ f dT[U mm ,n^n m ' a + (U mml -J mm ,)nZ l n^\ (35) 

m,m' ,<r ^ 

It describes the spin, orbital, energy and temperature dependent inter- 
actions of a particular magnetic 3(i-atom with the rest of the crystal and 
is used to compute the local Greens function matrix: 

G CT (r - T ') = -^J D[c,c + ]e- s ^c(r)c + (r') (36) 

[Z is the partition function) and the impurity self-energy Q~ l (uj n ) — 

G-Vn) = S*K) ■ 

The Weiss field function is required to obey the self-consistency con- 
dition (14), which can be specified for a given case as 

G CT K) = ]T[(*u, n + t i)i- n LDA (k) - s^K)]- 1 (37) 

k 

The local matrix X!^ c is the sum of two terms, the impurity self energy 
and a so-called "double counting " correction, E^c which is meant to 
subtract the average electron-electron interactions already included in 
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the LDA Hamiltonian. For metallic systems we proposed the general 
form of dc-correction: X!^ c (u) = S CT (u) — \Tr a Y, a (0). This is motivated 
by the fact that the static part of the correlation effects are already 
well described in the density functional theory Only the d-part of the 
self-energy is presented in our calculations, therefore £^ c = for s- 
and p- states as well as for non-diagonal d — s, p contributions. In 
order to describe the finite temperature ferromagnetism of transition 
metals we use the non spin-polarized LDA Hamiltonian H LDA (k) and 
accumulate all temperature-dependent spin-splittings in the self-energy 
matrix (o> n ). 

We used the impurity QMC scheme (see Appendix) for the solution 
of the multiband DMFT equations [45]. The Hirsch discrete Hubbard- 
Stratonovich transformations [43, 44] introduces (2M — 1)M auxiliary 
Ising fields , where M is the orbital degeneracy of the cZ-states and 
calculate G CT (r) by an exact integration of the fermion degrees of free- 
dom in the functional integral (Eq.(36)) [3]. In order to sample efficiently 
all the spin configurations in the multi-band QMC scheme, it is impor- 
tant to use "global" spin-flips: [S™^,] -> [-S~^ a '] in addition to the 
local moves of the auxiliary fields. The number of QMC sweeps was of 
the order of 10 5 . A parallel version of the DMFT program was used to 
sample the 45 Ising fields for 3d-orbitals. We used 256 k-points in the 
irreducible part of the BZ for the k integration. 10 to 20 DMFT iter- 
ations were sufficient to achieve convergence far from the Curie point. 
Due to the cubic symmetry of the bcc-Fe and fcc-Ni lattices the local 
Green function is diagonal in the basis of real spherical harmonics. The 
spectral functions for real frequencies were obtained from the QMC data 
by applying the maximum entropy method [46]. 

Our results for the local spectral function for iron and nickel are shown 
in Figs.l and 2, respectively. The SDF approach describes well all the 
qualitative features of the density of states (DOS), which is especially 
non-trivial for nickel. Our QMC results reproduce well the three main 
correlation effects on the one particle spectra below Tc [47, 48, 49]: the 
presence of a famous 6 eV satellite, the 30% narrowing of the occupied 
part of ci-band and the 50% decrease of exchange splittings compared 
to the LDA results. Note that the satellite in Ni has substantially more 
spin-up contributions in agreement with photoemission spectra [49]. The 
exchange splitting of the d-band depends very weakly on temperature 
from T=0.6 Tc to T=0.9 Tc. Correlation effects in Fe are less pro- 
nounced than in Ni, due to its large spin-splitting and the characteristic 
bcc-structural dip in the density of states for spin-down states near Fermi 
level, which reduces the density of states for particle hole excitations. 
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Now we discuss the applications of the SDF approach to the descrip- 
tion of the finite-temperature magnetic properties of iron and nickel. 
While density functional theory can in principle provide a rigorous de- 
scription of the thermodynamic properties, at present there is no accu- 
rate practical implemention available. As a result the finite-temperature 
properties of magnetic materials are estimated following a simple sug- 
gestion [30], whereby constrained DFT at T = is used to extract 
exchange constants for a classical Heisenberg model, which in turn is 
solved using approximation methods (e.g. RPA, mean field ) from clas- 
sical statistical mechanics of spin systems [30, 50, 51, 52]. The most 
recent implementation of this approach gives good values for the tran- 
sition temperature of iron but not of nickel [53]. While these localized 
spin models give, by construction, at high temperatures a Curie- Weiss 
like magnetic susceptibility, as observed experimentally in Fe and Ni, 
they encounter difficulties in predicting the correct values of the Curie 
constants [54]. 

The uniform spin susceptibility in the paramagnetic state, Xq=o = 
dM/dH, was extracted from the QMC simulations by measuring the in- 
duced magnetic moment in a small external magnetic field. It includes 
the polarization of the impurity Weiss field by the external field [3] . The 
dynamical mean field results account for the Curie- Weiss law which is ob- 
served experimentally in Fe and Ni. As the temperature increases above 
Tc, the atomic character of the system is partially restored resulting in 
an atomic like susceptibility with an effective moment: 

^2 

Xq=0 = 3(T -T C ) (38) 

The temperature dependence of the ordered magnetic moment below the 
Curie temperature and the inverse of the uniform susceptibility above 
the Curie point are plotted in Fig. 3 together with the corresponding 
experimental data for iron and nickel[55]. The LDA+DMFT calculations 
describes the magnetization curve and the slope of the high-temperature 
Curie- Weiss susceptibility remarkably well. The calculated values of 
high-temperature magnetic moments extracted from the uniform spin 
susceptibility are /x e // = 3.09 (1.50)//b for Fe (Ni), in good agreement 
with the experimental data fi e ff = 3.13 (1.62)/i£ for Fe (Ni)[55]. 



We have estimated the values of the Curie temperatures of Fe and Ni 
from the disappearance of spin polarization in the self-consistent solu- 
tion of DMFT problem and from the Curie- Weiss law in Eq.(38). Our 
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Figure 1. LDA+DMFT results for ferromagnetic iron (T = 0.8 Tc)- The partial 
densities of d-states (full lines) is compared with the corresponding LSDA results at 
zero temperature (dashed lines) for the spin-up (arrow-up) and spin-down (arrow- 
down) states. The insert shows the spin-spin autocorrelation function for T=1.2 Tc- 
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Figure 2. Same quantities as in Fig.l for ferromagnetic nickel (T = 0.9 Tc). The 
insert shows the spin-spin autocorrelation function for T=1.8 Tc- 
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Figure 3. Temperature dependence of ordered moment and the inverse ferromag- 
netic susceptibility for Fe (open square) and Ni (open circle) compared with experi- 
mental results for Fe (square) and Ni (circle) (from Ref.[33]). The calculated moments 
were normalized to the LDA ground state magnetization (2.2 /j,b for Fe and 0.6 fj,B 
for Ni). 

estimates Tq = 1900 (700) if are in reasonable agreement with experi- 
mental values of 1043 (631) K for Fe (Ni) respectively [5 5], considering the 
single-site nature of the DMFT approach, which is not able to capture 
the reduction of Tq due to long wavelength spin- waves. These effects 
are governed by the spin- wave stiffness. Since the ratio of the spin- wave 
stiffness (D) to Tc, Tc/a 2 D is nearly a factor of 3 larger for Fe than for 
Ni[55] (a is the lattice spacing), we expect the DMFT Tc to be much 
higher than the observed Curie temperature in Fe than in Ni. Note that 
this is a consequence of the long-range oscilating character of exchange 
interactions in iron compared to short-range ferromagnetic exchange in- 
teractions in nickel [53]. Quantitative calculations demonstrating the 
sizeable reduction of Tq due to spin waves in Fe in the framework of a 
Heisenberg model were performed in Ref [53] . 

Within dynamical mean field theory one can also compute the local 
spin susceptibility defined by 




(39) 



o 
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where g s = 2 is the gyromagnetic ratio and S = \ ^2 mcra i v acr > >c ma i is 
single-site spin operator and a = (cr x ,a y ,a z ) are Pauli matrices. It dif- 
fers from the q = susceptibility by the absence of spin polarization in 
the Weiss field of the impurity model. Eq.(39) cannot be probed directly 
in experiments but it is easily computed in DMFT-QMC. Its behavior 
as function of temperature gives a very intuitive picture of the degree 
of correlations in the system. In a weakly correlated system we expect 
Eq.(39) to be nearly temperature independent, while in a strongly corre- 
lated system we expect a leading Curie- Weiss behavior at high tempera- 
tures Xiocal = fJ-foc/ fiT + const) where ui oc is an effective local magnetic 
moment. In the Heisenberg model with spin S, fsf oc = S(S+l)gg and for 
well-defined local magnetic moments (e.g., for rare earth magnets) this 
quantity should be temperature independent. For the itinerant electron 
magnets, \i\ oc is temperature-dependent, due to a variety of compet- 
ing many-body effects such as Kondo screening, the induction of local 
magnetic moment by temperature [35] and thermal fluctuations which 
disorders the moments [56]. All these effects are included in the DMFT 
calculations. The r-dependence of the correlation function (S(r)S(O)) 
results in the temperature dependence of m oc and is displayed in the 
inserts on the Figs. 1,2. Iron can be considered as a magnet with very 
well-defined local moments above Tq (the r-dependence of the correla- 
tion function is relatively weak), whereas nickel is more of an itinerant 
electron magnet (stronger r -dependence of the local spin-spin autocor- 
relation function). 

The comparison of the values of the local and the q = susceptibil- 
ity gives a crude measure of the degree of short-range order which is 
present above Tc- As expected, the moments extracted from the local 
susceptibility Eq.(39) are a bit smaller ( 2.8 fiB for iron and 1.3 fiB 
for nickel) than those extracted from the uniform mangetic susceptibil- 
ity. This reflects the small degree of the short-range correlations which 
remain well above Tc [58]. The high-temperature LDA+DMFT clearly 
show the presence of a local- moment above Tc- This moment, is cor- 
related with the presence of high energy features (of the order of the 
Coulomb energies) in the photomeission. This is also true below Tc, 
where the spin dependence of the spectra is more pronounced for the 
satellite rigion in nickel than for that of the quasiparticle bands near 
the Fermi level (Fig. 2). This can explain the apparent discrepancies 
between different experimental determinations of the high-temperature 
magnetic splittings [57, 59, 60] as being the results of probing differ- 
ent energy regions. The resonant photoemission experiments [59] reflect 
the presence of local-moment polarization in the high-energy spectrum 
above the Curie temperature in nickel, while the low-energy ARPES 



SDF approach to electronic correlations 



19 



investigations [60] results in non-magnetic bands near the Fermi level. 
This is exactly the DMFT view on the electronic structure of transition 
metals above Tq- Fluctuating moments and atomic-like configurations 
are large at short times, which results in correlation effects in the high- 
energy spectra such as spin-multiplet splittings. The moment is reduced 
at longer time scales, corresponding to a more band-like, less correlated 
electronic structure near the Fermi level. 

7. Approximate solution for the self-energy: 



The QMC method described above is probably the most accurate 
way of solving the effective impurity problem in the DMFT. However, 
it is rather cumbersome and expensive computationally; besides that, it 
deals with the "truncated" two- indices interaction matrix (see Eq.(35)) 
instead of the complete four-indices one. Therefore a scheme has been 
proposed in Ref. [6] based on a multiband spin-polarized generalization 
of the "fluctuating exchange" (FLEX) approximation by Bickers and 
Scalapino [61]. The original formulation of the FLEX approximation 
treats both particle-hole (PH) and particle-particle (PP) channels on an 
equal footing. But their roles in magnetism are completely different. 
The interaction of electrons with spin fluctuations in PH channel leads 
to the most relevant correlation effects [35] whereas PP processes are 
important for the renormalizations of the effective interactions in spirit 
of the T-matrix approach ("ladder aprroximation" ) by Galitskii [62] and 
Kanamori [63]. Therefore we used in Ref. [6] a "two-step" procedure 
when, at first, the bare matrix vertex is replaced by a T-matrix, and, 
secondly, PH channel processes with this effective interaction are taken 
into account explicitly. This approximation has high enough accuracy 
both for the Hubbard model and for real systems with moderate cor- 
relations U < W/2 where U is the Hubbard on-sire repulsion snergy 
and W is the bandwidth (see [6] and Refs therein). However, specific 
form of the approximation used in [6] can be improved further by taking 
into account the spin-dependence of the T-matrix. Here we present the 
formulation of this T-matrix-FLEX approximation. 

Let us start, as in Refs. [5, 6], with the general many-body Hamilto- 
nian for a crystal in the LDA+U scheme [64]: 



TM-FLEX method 



H = H t + Hu 
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{\i}<ra' 



(AiA 2 |v| AiAa) 4 1CT C A 2 a' C A^'C A ' lCT , 



(40) 



where A = im are the site number (i) and orbital (m) quantum numbers, 
a =|, J, is the spin projection, c + , c are the Fermion creation and annihi- 
lation operators, Ht is the effective single-particle Hamiltonian from the 
LDA, corrected for the double-counting of average interactions among 
correlated electrons as it was described above, and the Coulomb matrix 
elements are defined in the standard way 

(12 |v| 34) = J drdr^l(r^*(r')v (r - r') ^(r)^^'), (41) 

where we define for briefness Ai = 1 etc. Following Ref. [62] we take 
into account the ladder (T-matrix) renormalization of the effective ap- 
proximation: 



13 



T aa ' (in) 24) = (13 \v\ 24)- ^EE < 13 M 57 ) * 



P w 5678 

Gl 6 (iu) G^' 8 (in - iu) (68 It™' (in) 



24 



(42) 



Further we rewrite the perturbation theory in terms of this effective 
interaction matrix. 

At first, we take into account the "Hartree" and "Fock" diagrams with 
the replacement of the bare interaction by the T-matrix 



Eg? (H 



= -^EE( 14 I TCTCT (^)I 32 ) G 34(^-^) (43) 



\ E E ( 13 \ T<Ta ' H 24 > G « w - ™ 



U 34<t' 
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CI 34 



Note that T,( TH ^+ Y,( TF ^ contains exactly all the second-order contri- 
butions. Now we have to consider the contribution of particle- hole ex- 
citations to sigma. Similar to [6] we will replace in the corresponding 
diagrams the bare interaction by the static limit of the T-matrix. How- 
ever, we improve the approximation [6] by taking into account its spin 
dependence. When considering the particle-hole channel we replace in 
the Hamiltonian (40) v -> T CT<T ' which is the solution of Eq.(42) at n = 0. 
Eq. (43) is exact in the limit of low electron (or hole) density which is 
important for the criterion of magnetism e.g. in the case of nickel (with 
almost completely filled d- band). 

Now we rewrite the effective Hamiltonian (40) with the replacement 
(12 |u| 34) by (l2 T aa '(n = 0) 34) in H v . To consider the correlation 
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effects due to PH channel we have to separate density (d) and magnetic 
(m) channels as in [61] 



m 12 



-/= (c| T c 2T - c+ c 2i ) 



m 



12 



m 12 = 



C U C 2T 



(44) 



Then the interaction Hamiltonian can be rewritten in the following ma- 
trix form 

H v = ^Tr [D + * VH * D + m+ * v£ * mT + m~ * * m+) (45) 

where * means the matrix multiplication with respect to the pairs of 
orbital indices, e.g. 



34 

the supervector D defined as 

D = (* m »),^_(^), 
and the effective interactions have the following form: 

(v£) 

V m /ll',22' 
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(46) 
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To calculate the PH contribution to the electron self-energy we first 
have to write the expressions for the generalized susceptibilities, both 
transverse x 1 ' an d longitudinal x'L The corresponding expressions are 
the same as in [6] but with another definition of the interaction vertices. 
One has 



1 *T n (iu), (47) 
where 

n( 34 (r) = -G^(r)G^(-T) (48) 

is an "empty loop" susceptibility and T(iu) is its Fourier transform. The 
corresponding longitudinal susceptibility matrix has a more complicated 
form: 



IUJ 



x"M = [l + y» * x8(tw)] *x W), 

and the matrix of bare longitudinal susceptibility: 



(49) 



Xo- 2 



i f r^T + r u r TT - r u 

pTT _ pU pTT _|_ pU 



(50) 



in the dd-, dm -, m°d-, and m°m - channels (d, m° = 1, 2 in the super- 
matrix indices). An important feature of these equations is the coupling 
of longitudinal magnetic fluctuations and of density fluctuations. It is 
absent in one-band Hubbard model due to the absense of the interaction 
of electrons with parallel spins. For this case Eqs. (47,49) coinsides with 
the well-known result [65]. 

Now we can write the particle-hole contribution to the self-energy. 
According to [6] one has 



Sffi (r) = E W is% (t) G° u (t) , 

34,(7' 



(51) 



with P-H fluctuation potential matrix: 
W aa ' {iuj) -- 



W n (iuj) W 1 - (iu) 
W ± {iuj) W 11 (iio) 



(52) 



where the spin-dependent effective potentials are defined as 
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where x">Xo differ from x">Xo by the replacement of 44> in 
Eq.(50). We have subtracted the seconf-order contributions since they 
have already been taken into account in Eq.(43). 
Our complete expression for the self energy is 



(53) 



This expression takes into account accurately spin-polaron effects be- 
cause of the interaction with magnetic fluctuations [6], the energy de- 
pendence of T-matrix which is important for describing the satellite 
effects in Ni [36], contains the exact second-order terms in v and is rig- 
orous (because of the first term) for almost filled or almost empty bands. 
In spirit of the DMFT approach we have to use Qq instead of G in all 
the expressions when calculating the self-energy on a separated central 
site. It should be noted that this TM-FLEX scheme is not conserved (or 
"^-derivable" ) therefore one need to inforce the Luttinger theorem by 
introducing the (jlq for the bath Green function as in iterative preturba- 
tion theory [66] 

We have started from the spin-polarized LSDA band structure of fer- 
romagnetic nickel within the TB-LMTO method [41] in the minimal 
s,p,d basis set and used numerical orthogonalization to find the Ht part 
of our starting Hamiltonian. We take into accounts the Coulomb inter- 
actions only between d-states. Semiempirical analysis of the appropriate 
interaction value gives U ~ 2 — 4 eV. The difficulties with choosing the 
correct value of U are connected with complicated screening problems, 
definitions of orthogonal orbitals in the crystal, and contributions of the 
intersite interactions. In the quasiatomic (spherical) approximation the 
full [/-matrix for the d— shell is determined by the three parameters U, J 
and 5 J or equivalently by effective Slater integrals F , F 2 and F 4 [64]. 
For example, U = F°, J = (F 2 + F 4 )/14 and we use the simplest way of 
estimating 5 J or F 4 keeping the ratio F 2 /F 4 equal to its atomic value 
0.625 [64]. Note that the value of intra-atomic (Hund) exchange inter- 
action J is not sensitive to the screening and approximately equals 0.9 
eV in different estimations [64]. 
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Figure 4- Spin- up (full lines) and spin-down (dased lines) density of d-states for fer- 
romagnetic nickel in the LSDA and the LDA+SPTF (LSDA+SPTF) calculations for 
different average Coulomb interaction U with J = 1 eV and temperature T=200 K. 
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The resulting DOS for Ni (Fig. 4) shows that spin-polarized TM- 
FLEX calculations approximately reproduce the satellite structure and 
reduction of the band width in satisfactory agreement with exact QMC- 
result(Fig. 2) 

8. Cluster DMFT appoach: antiferromagnetism 
and superconductivity 

As it was stressed above the SDF approach does not necessarily con- 
nect with the standard DMFT scheme; one can use more general choise 
of fermion coarse-grained variables, e.g., considering the case of an effec- 
tive cluster instead of an effective impurity. It is especially important for 
the problems where intersite correlations are involved from the begin- 
ning such as d-wave superconducting pairing [22, 23] or charge ordering 
[67]. Here we consider, following Ref. [23], one of the cluster general- 
izations of the DMFT, and its application to the problem of magnetism 
and high-temperature superconductivity (HTSC) of copper-oxide com- 
pounds. 

The microscopic theory of high-temperature superconducting cuprates 
is still far from a final understanding [68, 69, 70]. One of the most impor- 
tant recent experimental achievements was the discovery of the pseudo- 
gap (PG) phenomenon above the superconducting transition tempera- 
ture [71] and existence of a sharp 41-meV resonance below T c related to 
some collective antiferromagnetic excitations [72]. Thus, an interplay of 
an antiferromagnetism (AFM) and <i-wave superconductivity (d-SC) in 
cuprates could be a natural way of discussing different HTSC phenom- 
ena. This require a quantitative electronic structure theory including 
two different type of the order parameters: AFM and d-SC. Within such 
approach one can in principle analyze the phase diagram of HTSC com- 
pounds and resolve the long-standing problem of competition between 
antiferromagnetism and d-wave superconductivity in cuprates [73, 74]. 

A standard theoretical tool for cuprates electronic structure consists 
of the two-dimensional Hubbard model [68] . We start with the extended- 
hopping Hubbard model on the square lattice: 

H = J2 tij4a c jv + J2 U i n i\ n H 

ij i 

where Uj is an effective hopping and Ui local Coulomb interactions. We 
have chosen the nearest-neighbor hopping t = 0.25 eV and the next 
nearest hopping t' /t = —0.15 for the model of La2- x Sr x Cu04[75]. The 
total band width is W=2 eV and all Coulomb parameters set to be 
£7=1.2 eV (U/W = 0.6). Let us introduce the "super-site" as an 2x2 
square plaquet. The numeration of the atoms in the super-site is shown 
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in the Fig.5. It is useful to introduce the superspinor Cf = {cf a } where 
a = 0, 1, 2, 3 (the spin indices are not shown). Taking into account the 
spin degrees of freedom, this is the 8-component superspinor creation 
operator. Then the crystal Green function for the Hubbard model can 
be rewritten as 



G (k,iu) = [iuj + fi - h (k) - £(zu;)]" 



(54) 



where h (k,zu>) is the effective hopping supermatrix with self-energy cor- 
rections. For simplicity we will write all the formulas in the nearest- 
neighbor approximations: 



h(k) 



( 

tK: 







tK+ \ 
tK+ 







tK~ tK x 
tK~ tK+ J 



(55) 



where K^, , 



1 + exp (±ik x ( y -)Ctj , a is the lattice constant, and each ele- 
ment is a 2x2 matrix in spin space. Within the cluster-DMFT approach 
we introduce the intra-atomic self-energy So and the inter-atomic self- 
energies S x , Tj y , and both functions are of intra-site nature in the sense 
of our super-site: 



S (iu>) 



So 
S* 


s* 



So 






Sj, 
So 

S. r 



s, \ 



s* 

So / 



(56) 



For the small 2x2 cluster it is usefull to introduce the translation- 
ally invariant (k-dependent) self-energy and rewrite h (k) + S(iw) — > 
h (\t,iuj) where 



/ S C 



h Qc,iu) = 



t*K~ 



t K + 





t K + 



t K + 


So 



\ 



(57) 



The effective Hamiltonian is defined through the renormalized energy 
dependent hoppings: t x = t + S x , t y = t + T, y . The functions So (iuj) , 
£ x (ico) , Tj y (iio) are found self-consistently within the cluster DMFT 
scheme and for the d-wme superconduction state S x ^ T, y . It is straight- 
forward to generalize this scheme for a next-nearest neighbor hopping as 
well as the long-range Green function and the self-energy. In this case 
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Figure 5. (a) Schematic representation of antiferromagnetic d-wave 2x2 periodically 
repeated clust er; (b) generic phase diagram of HTSC materials; (c) The calculated 
values of two order parameters: local magnetic moment M and d-SC equal time Green 
function F 01 (t = 0) = F(0) for different hole doping (x) at the inverse temperature 
/3 = 60 eV _1 (T = 190K). 



we can renormalized also the second-nearest hopping: t 



xy 



t' + Ti xy for 



the 2x2 cluster, where Y, xy (or E ) is the non-local self-energy in the 
xy direction. 

In the cluster version of the DMFT scheme one can write the matrix 
equation for the bath Green function matrix Q which describe an effective 
interaction with the rest of crystal: 

G~ l (iu) = G" 1 (ioj) + E (iu) , 

where the local cluster Green function matrix is equal to G a /3 (iu) = 
S G a fi (k,iu) , and the summation is run over the Brillouin zone of 

k 

the square lattice. If instead of Eq.(57) we use Eq.(55), it would cor- 
responds to a free cluster DMFT scheme or the simplest case of the 
so-called cellular-DMFT. Note that we use a translationally invariant 
self-energy obtained from the cluster DMFT scheme (Eq.(57)) or the 
simplest version of so-called cellular DMFT approach [24]. We believe 
that for a given choise of the small (2x2) cluster, the renormalization 
of the hopping parameter by E x , T, y terms is physically essential. The 
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present "matrix" form of a cluster DMFT with the self-energy which 
is not periodic inside the cluster allow us to study a multicomponent 
ordered state. Unfortunately, in contrast with the so-called Dynamical 
Cluster Approximation (DCA) [22] or cellular-DMFT approach [24] (see 
above, Section 1), we were unable to prove the casuality of this approach 
for arbitrary band structure and interaction parameters. However, the 
casuality of the Green function for the realistic choise of the parameters 
has been checked numerically. 

In this case we have the standard DMFT problem with four "orbital" 
states per super-site. It has been solved by the multi-orbital QMC tech- 
nique described above (Section 3). We used the generalized Nambu 
technique [76] to analyze the coexistence of the magnetic ordering and 
superconductivity. Let us introduce the superspinor 

and the anomalous averages describing the (collinear) antiferromagnetism 
c V\ c ji) an d the superconductivity Ajj = (c^c^). 
The generalization of the Hirsch-Fye QMC-algorithm [44] for the su- 
perconducting problem [77] have been used. In the 4-spinor case a dis- 
crete Hubbard-Stratonovich transformation has the following form: 



exp[-ArL^nj T ni| H ^( n «t + n «l)l = 

\ exp[AiO-(V>H^ii - ~ $j#3i + tpiiipii)} (58) 

Z a=±l 

where Aj = ^arccosh[exp(iArC/j)]. 

Since we take into account only the singlet pairing, we obtain the 
following nonzero elements of the d-SC energy gap parameters: A = 
Aoi = — A12 = A23 = — A30. One can chose Ajj to be real and therefore 
symmetric: Ajj = Ajj. Separating normal and anomalous parts of the 
Green function we have 

^-{^&) <■»> 

where G (k,r, r') = - (T T C k (r) C+ (r')) and F (k,r, r') = - <T T C k (r) C_ k (r')> 
are the matrices in spin and "orbital" space. It is convenient to expand 
the anomalous Green function in Pauli matrices F = (F° + Fa) ia y and 
use the symmetry properties [78]: 

F°(k,r,r') = F°(-k,r',r) ( 60 ) 

F(k, t,t') = -F(-ky,r). 
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then a 4 x 4 spinor formalism is reduced to 2 x 2 one in the collinear anti- 
ferromagnetic case with the d-wave superconductivity with the following 
spin-matrix form of the local Green function for the super-site: 




and the QMC formalism for the antiferromagnetic superconducting state 
is equivalent to the previous non-magnetic one [77] . Using the discretiza- 
tion of [0,(3] interval with L-time slices: At = (3/L {(3 = 1/T is an in- 
verse temperature) the G a - and F- Green functions become the matrices 
of 2NL dimension, where N is the number of atoms in the cluster. Af- 
ter Fourier transform to the Matsubara frequencies the Green function 
matrix has the following form: 

G ^)-y F{iu) -G\{iu) ) (62) 
In superconducting states the self-energy defined as [3] : 

G M-G ( W )=^ ^ M _ E j (jw) j, (63) 
and the inverse crystal Green function matrix is equal to: 

\ s (k,zw) iuj - n + h* (k,zw) J v ' 

where s (k,iu;) is the translationally invariant anomalous part of the self- 
energy S(iu) similar to Eq.(55). 

The two-component order parameters state which includes Neel anti- 
ferromagnetism and d-wave superconductivity (Fig. 5a) lowered the sym- 
metry of the effective cluster-DMFT problem. A self-consistent DMFT 
cluster problem with AFM and d-SC general order parameters have been 
solved within the QMC scheme for 8x8 matrix Green function with L=64 
time slices. The resulting two order parameters for (3 = 60 eV _1 (T=190 
K) and if = presented in Fig.5c together with the generic HTSC phase 
diagram (Fig. 5b) as function of the hole doping. In this case the or- 
dered magnetic moment is directly related with imaginary-time Green 
function G a (t) : M = G^° (0) - Gf (0) and for the d-SC order parame- 
ter we chose a positive value of superconducting imaginary time Green 
function F 01 (0). It is important, that we find no serious sign-problem 
for all QMC calculations with various doping level, probably due to "sta- 
bilized" antiferromagnetic dynamical mean fields acting on the atoms in 



30 




x (eV" ) 



Figure 6. Imaginary time normal (G CT ) and superconducting (F) Green functions 
for the 2x2 cluster DMFT solution with second-nearest neighbor hopping and inverse 
temperature (3 = 50 eV" 1 (T=230 K). 



our 2x2 cluster. Note that the AFM cluster-DMFT solution exists for a 
much higher doping concentration than experimental AFM ordered state 
and describes a dynamical mean-field version of AFM-spin fluctuations 
related to pseudogap phenomena (the PG-region on Fig. 5b). The max- 
imum of d-SC order parameter corresponds to a doping level of about 
15% in agreement with the generic HTSC phase diagram. The d-SC or- 
der parameter is zero close to the undoped region (x=0), due to presence 
of a large AFM-gap. When the AFM-gap is closed (x ~ 5%) the d-SC 
states develope but for x > 20% the d-gap decrease again since AFM 
spin-fluctuations around (vr,7r) point disappear [69]. The precise char- 
acteristic of the phase diagram including the interactions between the 
AFM and d-SC order parameters demands an extensive cluster-DMFT 
calculations for different temperatures and doping. 

We would like to note that the existing of d-SC cluster-DMFT solution 
for such high temperatures does not necessary means that the supercon- 
ducting transition temperature is larger then 190K in our model. A 
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crude estimation shows that the d-SC solution disappears at T=300 K 
for £=0.15 and the AFM solutions for x=0 become unstable at the tem- 
perature just above 1000 K. This could be the sign of a "local" AFM 
solution and a local d-wave solution, like local moments in magnetic 
systems [6]. Due to the multiscale nature of the problem under con- 
siderations, essentially different energies connected with local moment 
formation, long-range magnetic order, local d-wave pairs within the 2x2 
plaquet, and finally coherent superconductivity, it is difficult to distin- 
guish a real long-range ordering from slow dynamical fluctuations in our 
QMC simulations. We plane to separate these energy scales analyti- 
cally and estimate superconductiong transition temperature in a future 
publication. 

The role of next-nearest hopping is to lower of the van-Hove singu- 
larity [75] which increases the density of state at the Fermi level for the 
hole-doped case and favores the d-SC solution for a moderate correla- 
tion strength. There is also a change in the spin-fluctuation spectrum 
related with the broadening of AFM-peak near the (tt, tt) point due to 
formation of so-called extended van-Hove singularities with increasing 
of the t'. We show one of the AFM-dSC solution on the Fig. 6 with 
the next nearest-neighbor hopping for the 10% doping level and (3 = 50 
cV _1 . The resulting local magnetic moment is M=0.28 and the d-SC 
order parameter F(0)=0.036. One can see that the superconduction or- 
der parameter is really of the d x 2_ y 2 symmetry since diagonal elements 
(F 00 ) as well as the next nearest-neighbors elements (F 02 ) are all equal 
to zero and only the nearest neighbor superconducting Green functions 
(F 01 ) are non-zero and change the sign for F x and F y components. The 
normal local Green function (G 00 ) (plotted for the spin-up atom in the 
Fig.6) as well as (G 02 ) are spin-split, while the nearest-neighbor Green 
function (G 01 ) has no spin-splitting due to AFM spin symmetry (see 
Fig.5). The absence of magnetic polarization in the non-diagonal G- 
function along the x(y)-directions suppress the magnetic pair-breaking 
and makes the AFM-dSC coexistence possible. 

9. Summary 

The spectral density functional (SDF) approach allows us to study 
correlation effects in solids based on realistic electronic structure calcula- 
tions. Among all possible applications we have chosen the magnetism of 
transition metals and high-temperature superconductivity. From these 
two examples one can see already all the main advantages of the new 
approach in comparison with standard density functional theory. First, 
we can describe now the spectral density transfer phenomena (e.g., the 
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formation of 6 eV satellite in Ni), the quasiparticle damping and other 
effects connected with the frequency dependence of the self-energy; they 
are absent completely not only in the DF approach but also in the 
Hartree-Fock, LDA+U, or self-interaction corrections approximation (see, 
e.g., [64]). Second, we can describe adequately the contribution of the 
Bose degrees of freedom (e.g., spin fluctuations) to the electronic struc- 
ture and thermodynamic properties. In the DF-based calculations the 
temperature is really taken into account only via the thermal expan- 
sion and the Fermi distribution function [79]. It was the main reason 
for the failure of the standard band theory for the description of finite- 
temperature effects in magnetic metals. We show that the SDF gives a 
satisfactory solution of this problem. 

Most of real applications of the SDF approach are connected with the 
single-site dynamical mean field theory. At the same time, for a number 
of problems this can be insufficient and some generalizations to take into 
account the non-local effects are necessary, for example, cluster ones. 
The first attempts of such generalizations are already leading to some 
interesting results in the theory of high-T c superconductivity [22, 23] 
but a lot of additional work is required in the area of cluster dynamical 
mean field theories to reach the level of understanding that was reached 
in single site DMFT. 
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Appendix: Multi-orbital QMC scheme 

The multi-orbital DMFT problem and general cluster DMFT scheme can be re- 
duced to the general impurity action: 

Hi Hi 1 r-li 

s =- dr dr ' y2 c t( r )Si 3 ( T - T ') C A T ') + 7; / dTy^ni(T)Uijnj(T) 
Jo Jo id l Jo 

where i = {m, a} - orbital (site) and spin. Without spin-orbital coupling we have: 

The auxiliary fields Green-function QMC use the discrete Hubbard-Stratanovich 
transformation inroduced by Hirsch[43] 

explArUijlmnj - ^(ni + nj)]} = ^ ^ cxp{A i3 S' ij (n l - nj)} 

s«=±i 

where Sij (r) are the auxiliary Ising fields for each pair of orbitals and time slice with 
the strength: 

At 

Xij = arccosh[exp(-^-(7ij)] 

Using Hirsch transofrmation one can integrated out fermionic fields in the path in- 
tegrals] and resulting partition function and Green function matrix have the following 
form: 



r 

d = \^WJlY. GiS^detlG- 1 ^)} 

Si 3 (r) 

where JV/ is the number of Ising fields, L is the number of time slices, and G(Sij) is 
the Green function in the auxiliary Ising fields: 

A, = (e v '-l) 

here we introduce the generalized Pauli matrix: 

r + M < j 

> 3 

For efficient calculation of the Green function in arbitrary configuration of Ising 
fields Gij (S) we use the following Dyson equation [44] : 

v'-v 



G' = [1 + (1 -G)(e v - v - l)]"^ 
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The QMC important sampling scheme allowed us to integrate over the Ising fields 
with the abs(det[G _1 (£;./)] ) as a stochastic weight[44, 3]. For a single spin-flip Sij, 
the determinant ratio is calculated as following: 

det [G] I det [G'\ = R, R 3 - R l3 

Rr = l + [l-Gii(T,T)]Ai(T) 

Rj = l + [l-G^(r,r)]A,(r) 

Rij = G ij {T y T)A j (T)G ji {T,T)A i {T) 

and the Green function matrix updated in the standard maner[44, 3]: 



G' ilJ2 {T 1 ,T 2 ) = G tlj2 (T 1 ,T 2 ) + [G ni (n,-r) - 8 tl iS TUT ]Ai(r)/R i (T)G t]2 (T, r 2 ) 
G^(n,r 2 ) = G^, 2 (n,r 2 ) + [G^,(r 1; r) - 5 il ^ T1 , T ]A,(r)/Ji 3 .(r)G;, 2 (r,r 2 ) 
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